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Abstract 

We study the e ± a^ + El£ iss and l+l~ + E^ iss signatures (I = e, n) 
for different values of tan (5 in the mSUGRA model. With tan f3 rising, 
we observe a characteristic change in the shape of dilepton mass spec- 
tra in l + l~ + EJp iss versus e fjfi final states reflecting the decrease of 
X2 — * l + l~X± branching ratio. We also study the non mSUGRA mod- 
ifications of the CMS test point LM1 with arbitrary relations among 
gaugino and higgsino masses. For such modifications of the mSUGRA 
test point LM1 the number of lepton events depends rather strongly 
on the relations among gaugino and higgsino masses and in some mod- 
ifications of the test point LM1 the signatures with leptons and E™ lss 
do not lead to the SUSY discovery and the single SUSY discovery sig- 
nature remains the signature with n > 2 jets + EJp %ss + no leptons. 
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1 Introduction 



One of the goals of the Large Hadron Collider (LHC) pQ is the discovery 
of supersymmetry (SUSY). The squark and gluino decays produce missing 
transverse energy E™ tss from lightest stable superparticle (LSP) plus mul- 
tiple jets and isolated leptons pp. One of the most interesting and widely 
discussed signatures for SUSY discovery at the LHC is the signature with 
two opposite charge and the same flavour leptons [2]: l + l~ + E™ lss . The 
main reason of such interest is that neutralino decays into leptons and LSP 
X2 —> l + l~Xi contribute to this signature and the distribution of the l + l~ 
invariant mass rrii nv (l + l~) has the edge structure [3] that allows to determine 
some combination of the SUSY masses. 

The signature e ± / u T + E™ lss can be realized when X2 decays into r pair 
which is only relevant at large tan/3 [I). Also as it is shown in Ref. [3] at 
the level of CMSJET [5] simulation the use of e^y^ + EJp tss signature for 
large tan f3 allows to obtain nontrivial information on parameters of the decay 
X2 tt — > ttXi — > e y^y^yvvv. O n the other hand, the e ± /i T + E™ tss 
(with an arbitrary number of jets) can be used for the detection of lepton 
flavour violation in slepton decays [6], [7] at the LHC. 

In the Minimal Supersymmetric Model (MSSM) [8J supersymmetry is 
broken at some high scale M by generic soft terms, so in general all soft 
SUSY breaking terms are arbitrary which complicates the analysis and spoils 
the predictive power of the theory. In the Minimal Supergravity Model 
(mSUGRA) [8] the universality of the different soft parameters at the Grand 
Unified Theory (GUT) scale M GUT « 2 • 10 16 GeV is postulated. Namely, 
all the spin zero particle masses (squarks, sleptons, higgses) are postulated 
to be equal to the universal value tuq at the GUT scale. All gaugino particle 
masses are postulated to be equal to the universal value m.1/2 at GUT scale. 
Also the coefficients in front of quadratic and cubic SUSY soft breaking terms 
are postulated to be equal. The renormalization group equations are used 
to relate GUT and electroweak scales. The equations for the determination 
of a nontrivial minimum of the electroweak potential are used to decrease 
the number of the unknown parameters by two. So the mSUGRA model 
depends on five unknown parameters. At present, the more or less standard 
choice of free parameters in the mSUGRA model includes mo, m^, tan/3, A 
and sign(fi) [8]. All sparticle masses depend on these parameters. 

In Ref. [2] the possibility to detect SUSY and lepton flavour violation using 
the e^n^ + EJp tss signature at the LHC for the Compact Muon Solenoid 



1 



(CMS) detector at the level of full detector simulation was studied. Note 
also Ref.fTU] where the SUSY mSUGRA CMS detector discovery potential 
was investigated for the l + l~ + E™ lss signature at the level of full detector 
simulation. 

In this paper we study the e ± fi T + E^ tss and + E™ lss signatures 
(I = e, n) for different values of tan /3 in the mSUGRA model. With tan j3 ris- 
ing, we observe a characteristic change in the shape of dilepton mass spectra 
in +E™ lss versus e ± fi T final states reflecting the decrease of X2 ~> l + l~Xi 
branching ratio. We also study the non mSUGRA modifications of the CMS 
test point LM1 with arbitrary relations among gaugino and higgsino masses. 
For such modifications of the mSUGRA test point the number of lepton 
events depends rather strongly on the relations among gaugino and higgsino 
masses and in some modifications of test point LM1 the signatures with lep- 
tons and E™ lss do not lead to the SUSY discovery and the single SUSY dis- 
covery signature remains the signature with n > 2 jets + EJp iss + no leptons. 

The organization of the paper is the following. Section [2] describes some 
useful technical details of performed simulations. In Section [3] the back- 
grounds and cuts used to suppress the backgrounds are discussed. Section 
H] contains the results of numerical calculations concerning the possibility to 
detect SUSY and constrain the SUSY parameters using the e ± /i T + E™ lss 
and + E™ lss signatures. Section [5] contains the results of the simu- 

lations for different values of tan/?. In Section [6] we discuss the results of 
the simulation of the MSSM with nonuniversal gaugino and higgsino masses. 
Section [7] contains concluding remarks. 

2 Simulation details 

The coupling constants and cross sections in the leading order (LO) approx- 
imation for SUSY processes and backgrounds were calculated with ISAS- 
UGRA 7.69 [11J, PYTHIA 6.227 [12] and CompHEP 4.2pl rjLj. For the cal- 
culation of the next-to-leading order (NLO) corrections to the SUSY cross 
sections the PROSPINO [H] code was used. For considered signal events and 
backgrounds the NLO corrections are known and the values of NLO cross 
sections (or fc-factors) were used for normalization of the numerical results. 
We used the full simulation results of Ref . [9] for the estimation of the num- 
ber of background events. The CMS fast simulation code FAMOSAAD 
[T5] was used for the estimation of signal events. 
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The reconstructed electrons and muons were passed through packages 
defining lepton isolation criteria. For each electron and muon the following 
parameters were defined: 

• Trackl solation is a number of additional tracks with > 2 GeV/c 
inside a cone with R = \J Ar/ 2 + A$ 2 < 0.3 around the lepton. 

• C alol solation is a ratio of energy deposited in the calorimeters (elec- 
tromagnetic (ECAL) + hadronic (HCAL)) inside a cone with R = 0.13 
around given track to the energy deposited inside a cone with R = 0.3. 

• HEratio is defined as a ratio of energy deposited in the HCAL inside 
a cone with R = 0.13 to the energy deposited in the ECAL inside the 
same cone. 

• EPratio is a ratio of energy deposited in the ECAL inside a cone with 
R = 0.13 to the momentum of the reconstructed track. 



3 Signal selection and backgrounds 

The SUSY production pp — ► qq ,gg, qg with subsequent decays 

q^q'xi,2> (!) 
g -> qq'xtfr ( 2 ) 

Xt,2 - X?e + (/x>, (3) 

Xl,2 XifJ-~(e~)u, (4) 

lead to the event topologies e ± /i =F + and l + l~ + E'p tss . Note that 

in the MSSM with lepton flavour conservation neutralino decays into leptons 
X2 3 4 ~ * = e ; A 4 ) do not contribute into the e ± fi T + E™ lss signature. 

The main backgrounds which contribute to the e ± fi T events are: tt, WW, 
WZ, ZZ, Wt, Zbb, rf and Z+jet. It is found that tt is the largest background 
and it gives more than 50% contribution to the total background. In this 
paper we used the results of Ref. [5] for the estimation of the number of 
background events. 

In the analysis the events with the following isolation criteria for elec- 
trons were used: Trackl solation < 1.0, Calolsolation > 0.85, 0.85 < 
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EPratio < 2.0, HEratio < 0.25. The same criteria for muons were the 
following: Trackl solation < 1.0, Calolsolation > 0.50, EPratio < 0.20, 
HEratio > 0.70. These numbers were adjusted by studying electron and 
muon tracks in the process pp — > WW — ► 2/. 
The selection cuts are the following: 

• cut on leptons: p l ^ pt > p l ^ pt, °, \rj\ < 2.4, lepton isolation within AR < 
0.3 cone 

• cut on missing transverse energy: E™ tss > E^p lss ' . 
Here p l ^ pt '° and are corresponding thresholds. 

4 Use of the e ± fi T + i?™ ss signature for the 
SUSY detection 

In this section we remind the main results of Ref . [9] on SUSY detection using 
the e ± fr f + E™ tss signature. The possibility to detect SUSY using the CMS 
test points LM1 - LM9 [TB] chosen for the detailed study of SUSY detection 
at the CMS was investigated in Ref. |2] . This study was based on the counting 
the expected number of events for both the Standard Model (SM) and the 
mSUGRA model. The parameters of the CMS test points LM1 - LM9 are 
given in Tabled] 



Table 1: The parameters of the CMS test points. 



Point 


m (GeV) 


m 1 / 2 (GeV) 


tan (3 


sign(fi) 


A 


LM1 


60 


250 


10 


+ 





LM2 


185 


350 


35 


+ 





LM3 


330 


240 


20 


+ 





LM4 


210 


285 


10 


+ 





LM5 


230 


360 


10 


+ 





LM6 


85 


400 


10 


+ 





LM7 


3000 


230 


10 


+ 





LM8 


500 


300 


10 


+ 


-300 


LM9 


1450 


175 


50 


+ 
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For the point LM1 (the point LM1 coincides with the post-WMAP point 
B HZ!) it was found that the cuts with p l ^ pt > 20 GeV/c, E™ ss > 300 GeV 
are close to the optimal ones (the highest significance with the best sig- 
nal/background ratio). The results for the luminosity C = 10 fb~ are 
presented in Table H For other CMS SUSY test points LM2 - LM9 the 
results with the same cuts are presented in Table [3J The significances in 
this table are the following: S cl2 = 2(y / N s + N B — \fN B ) [TS] and S c l = 

y/2((N s + N B ) ln(l + N S /N B ) - N s ) [19]. 



Table 2: The expected number of events for backgrounds and for signal at 
the point LM1, £ = 10 fb' 1 , e ± ^ + EJp iss signature. 



Process 


2 isolated leptons, p l ^ pt > 20 GeV/c 


E™ iss > 300 GeV 


tt 


39679 


79 


WW 


4356 


4 


wz 


334 


2 


zz 


38 





wt 


3823 


2 


Zbb 


315 





Z+jet 


1082 


6 


DY2r 


7564 





SM background 


57191 


93 


LM1 Signal 


1054 


329 



It was found from the Tables [2J13] that for the point LM1 the significances 
are S c i 2 = 21.8 and S cL = 24.9 for the e ± /i T + E™ iss signature. 

The supersymmetry discovery potential for the mSUGRA model with 
tan/3 = 10, sign(fi) = + in the (m , m 1( / 2 ) plane (generalization of the 
point LM1) using the CMS fast simulation program FAMO 5*_1_4_0 [15] was 
also studied. The CMS discovery potential contours for C = 1, 10 and 30 fb^ 1 
for the signature e^fi 11 + E™ lss are shown in Fig{TJ 
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Table 3: The number of signal events and significances for the cuts with 
p l T Pt > 20 GeV/c and EJp iss > 300 GeV for £ = 10 fb~\ signature 
e^fi^ + EJp tss . The number of the SM background events Nb = 93 (see 
Table ED . 



Point 


iV events 


S C V2 


S C L 


LM1 


329 


21.8 


24.9 


LM2 


94 


8.1 


8.6 


LM3 


402 


25.2 


29.2 


LM4 


301 


20.4 


23.1 


LM5 


91 


7.8 


8.3 


LM6 


222 


16.2 


18.0 


LM7 


14 


1.4 


1.4 


LM8 


234 


16.9 


18.8 


LM9 


137 


11.0 


11.9 



5 The shape of the dilepton mass distribu- 
tion in the mSUGRA parameter space 

At the LHC neutralinos X2 are dominantly produced in the decay chain of 
gluino and squarks, for instance g — » qqx\ or Q ~> QX.2 an d Xj ~> X2Z . 
Within the mSUGRA model in the domain of relatively small mo (mo < 
0.6mi/ 2 ) the following leptonic decays of X2 are the most important: 

xi - x?J + r, (5) 

xi - y, (6) 

X°2 - hi, (7) 
X°2 - hr, (8) 

xl - x°ihn. (9) 

In mSUGRA for the case when X2 decays into lepton and slepton are 
kinematically allowed, the sleptons decay directly into LSP with Bt(Il 7 r 
lXi) ~ 1 and Br(?i — > txI) ~ 1. With the increase of tan/3 the fx mass 
is decreased and as a consequence the branching ratio Br(xX —> t\t) is 
increased whereas the branching ratio Br(xX ~~ * Ir,lI) is decreased, see Figj2j 
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The process —* U —> llXi h as the edge structure for the distribution 
of lepton-pair invariant mass m u and the edge mass r rn r ^{l + l~) is expressed 
by slepton mass mi and neutralino masses m^o ^ as follows [3] : 

m- m -o 

w(i + n) a = ^|o(i - -f )(i - -i) (io) 

z IIL~a llir 

xi i 

The corresponding Tn 1 ™ 1 J{T + T~) edge maximum due to x° decays to stau's 
has a maximum at: 

m 2 - m 2 ~o 

(C(^)f = m|o(l - -^)(1 - -£) (11) 



X2 



But the spectrum of the dilepton same flavour opposite sign channel pro- 
ceeding through t — ► \vv decays is not so pronounced as the spectrum from 
x\ — > 11 — > llXi decay due to the missing momemtum taken by four neutrinos 
from r decays. In particular the invariant mass is distributed in the lower 
mass region below the expected ditau kinematical point. Thus, a distinctive 
feature of the m inv {l + l ~) spectrum from decays to staus is distributed in 
the lower mass region below the expected ditau kinematical end point. With 
increase of tan (3, due to significant increase of x° decays into stau's and the 
corresponding decrease of decays into selectrons and smuons, a deterioration 
of the sharpeness of the dilepton edge takes place [I]. Another conse- 
quence of the change in the relative branching ratios of these decays is that 
the event rate difference between e + e~ + yU + /i~ and e + fi~ + fi + e~ channels 
decreases jl]. To illustrate the behaviour of the m inv (l + l~) and m inv (l + l ~) 
spectra at different tan (3 we investigated several mSUGRA points, namely: 

a. The generalization of the point LM1 (m = 60 GeV, m\/2 = 250 GeV, 
A = 0,sign(ji) = +,tan/3 = 10) with tan/5 = 15, 20, 25, 30, 35. 

b. The points with m = 100 GeV,m 1/2 = 300 GeV, A = 300 GeV, 
sign(fi) = + and tan/3 = 10, 15, 20, 25, 30, 35. 

c. The points with m = 300 GeV, = 300 GeV, A = 0, sign(fi) = + 
and tan/5 = 10, 15, 20, 25, 30, 35. 

d. The points with m = 2500 GeV, mi/2 = 250 GeV, A = 0, sign(fi) = + 
andtan/5 = 10, 15, 20, 25, 30, 35. 

For the points (a) and (b) sleptons are lighter than x\ that leads to the 
existence of well-defined edge structure in the m inv {l + l~) structure for not 
very large values of tan f3. For the points (c) and (d) the sleptons are heavier 
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than X2 an d the dominant decay mode is x\ ~ > % + Xi with the branching 
ratio close to 100%. 

For these mSUGRA points we have made calculations using FAMOS fast 
simulation program for the integrated luminosity C = 1 fb^ 1 . The results of 
our simulations are presented in FigsfJ] - [7J 

As one can see from Figsf2] - [7J with tan f3 rising there is a characteristic 
change in the shape of dilepton l + l~ mass spectra versus e^fj^ mass spectra 
resulting in the increase of the ratio of Ns(l + l~)/Ns(e ± fi T ) signal events. 
As it has been mentioned before this behaviour is due to the increase of 
Br{x2 ~^ an d decrease of Br(x2 ~ > H) with the rising of tan/3. 

6 The dileptons for the case of nonuniversal 
gaugino masses 

In this section we present the results of our simulations for the case of nonuni- 
versal gaugino masses. Namely, we consider the deformations of the test point 
LM1 violating gaugino mass universality assumption, see Table HI 



Table 4: Values of factors for masses increasing (unity if not indicated) in 
comparison with the masses of the point LM1. 



Point 


Mi 


M 2 


M 3 




1 




2 






2 




4 






3 




8 






4 




2 




2 


5 




4 




4 


6 




8 




8 


7 


2 


4 


5 


4 


8 


3 


4 


5 


4 


9 


2 


4 


5 


4 


10 


3 


4 


5 


4 


11 










12 


2 


2 







In Table H Mi, M 2 and M 3 are £7(1), SU(2) and SU(3) gaugino masses 



8 



correspondingly and /i is higgsino mass parameter. m 2 means masses of the 
first two squark and slepton generations and rriq means all squark and slepton 
masses. 

Points (1-3) differ from CMS mSUGRA test point LM1 by SU{2) gaug- 
ino mass M 2 . The dependence of rrii nv {l + l~) and mj rw (e ± /x =F ) invariant spec- 
tra for the point LM1 and the points (1-3) are shown in FigJSJ With the 
increase of SU(2) gaugino mass M 2 we see the sharp change in the dilepton 
invariant spectra. Points (4 - 6) differ from CMS mSUGRA test point LM1 
by both SU(2) gaugino mass and higgsino mass parameter \i. For the points 
(7 - 10) the difference with point LM1 also includes the increase in gluino 
and LSP masses. 

For the points (7 - 12) we have very sharp change in the number of dilepton 
events. For the luminosity £ = 1 fb^ 1 after the cuts with p l ^ pt > 20 GeV, 
fimiss > 3qq Q e y the number of events without isolated leptons, with single 
isolated lepton and with two isolated leptons are presented in Table [5j 

Table 5: The number of events with I — 0, I — 1, I — 2 for background and 
points LM1, 1-12^ 



Point N(l = 0) N(l = 1) N(l = 2) 



Back. 


20 


524 


652 


LM1 


2119 


1278 


166 


1 


1856 


1950 


430 


2 


1858 


1104 


187 


3 


1800 


671 


83 


4 


1910 


1876 


407 


5 


2002 


513 


34 


6 


1950 


449 


27 


7 


143 


35 


1 


8 


87 


21 





9 


95 


21 


1 


10 


56 


6 





11 


1587 


11 


1 


12 


1430 


10 


1 



As it follows from Table [5] for points (7 - 12) the perspective to discover 
SUSY using the signature with single lepton or dileptons plus E™ lss looks 



9 



hopeless. Moreover for the points (7 - 12) the use of more powerful from the 
SUSY discovery point of view signature no leptons + jets + E™ lss also 
looks very problematic. 

The invariant mass distributions of (solid line) and e^fi T (dotted 
line) lepton pairs at the extensions of the point LM1 are shown in FigsJB]- 
IHJ The distributions of leptons, 1 lepton, 2 leptons and > 3 leptons on 
Ejniss £ or moc iif] ca tions of the point LM1 are shown in Figs{TU]-[I3]for an 
integral luminosity C = 1 fb^ 1 . As one can see from Fig {10] for the points (7- 
10) there is sharp change in lepton spectra distributions on E™ tss compared 
to the test point LM1, namely, the number of lepton events for the points 
(7 - 10) are much smaller than for the test point LM1. The same situation 
takes place for the E™ tss spectrum of events without leptons, see Fig{TTJ 

7 Conclusion 

In this paper we studied the e^^ + E™ lss and l + l~ + E™ lss signatures (/ = 
e, /i) for different parameters tan (3 of the MSSM model. With tan f3 rising, 
we observed a characteristic change in the shape of dilepton mass spectra in 
+ EJp iss versus e^y^ final states reflecting the decrease of \2 ^ + ^ Xi 
branching ratio. We also studied some non mSUGRA modifications of the 
point LM1. We have found the drastic change in dilepton spectra with the 
increase of gaugino masses. For the points (7 - 12) the perspective to discover 
SUSY using the signatures with leptons and E™ lss looks hopeless. Even more 
powerful signature with no leptons + jets + E™ lss looks rather pessimistic. 
This work was supported by RFFI grant No 07-02-00256. 
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Figure 1: e + /i + e /i + discovery plot for tan/5 = 10, sign(fi) = +, 
A = 0. Two isolated leptons with p^ pt > 20 GeV/c and EJp iss > 300 GeV 
are selected. 
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Figure 2: The invariant mass distributions of l + l (solid line) and 
(dotted line) lepton pairs at the mSUGRA point m = 60 GeV, My 2 
250 GeV, [i > 0, A = with various tan/3 = 10, 15, 20, 25, 30, 35. 
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Figure 3: The invariant mass distributions of l + l (solid line) and e ± /i =F 
(dotted line) lepton pairs at the mSUGRA point m = 100 GeV,Mu 2 = 
300 GeV, [i > 0, A = 300 GeV with various tan/3 = 10, 15, 20, 25, 30, 35. 
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Figure 4: The invariant mass distributions of l + l (solid line) and e^ 1 
(dotted line) lepton pairs at the mSUGRA point m = 300 GeV, My 2 
300 GeV, fj. > 0,A = with various tan/5 = 10, 15, 20, 25, 30, 35. 
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Figure 5: The invariant mass distributions of l + l (solid line) and e 1 ^ 
(dotted line) lepton pairs at the mSUGRA point m = 2500 GeV, Mi/ 2 
250 GeV, n > 0,A = with various tan/5 = 10, 15, 20, 25, 30, 35. 
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Figure 6: The dependence of the ratio Ns(l + l~)/Ns(e ± fi T ) (Ns() is the 
number of events after the cuts p l ^ pt > 20 GeV, E™ tss > 300 GeV) on tan (5 
for different scenarios. Here the error is due to finite number of simulated 
events. 
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Figure 7: The dependence of the ratio Ns(l + l )/A r s(e ± /i T ) (Ns() is the 
number of events after the cuts p 1 ^ > 20 GeV, EJp iss > 300 GeV, \M inv - 
Mz\ > 15 GeV), on tan/? for different scenarios. Here the error is due to 
finite number of simulated events. 
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Figure 8: The invariant mass distributions of l + l (solid line) and e^[i 
(dotted line) lepton pairs at the extensions of the point LM1. 
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Figure 9: The invariant mass distributions of l + l (solid line) and e ± /i =F 
(dotted line) lepton pairs at the extensions of the point LM1. 
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Figure 11: The E™ tss distributions for the points LM1, 7 and 8 for the cases 
of no leptons and two leptons with p l ^ pt > 20 GeV. 
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Figure 12: The E™ tss distributions for the points LM1, 11 and 12. 
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Figure 13: The E^ tss distributions for points LM1, 11 and 12 for the cases 
of no leptons and two leptons with p l ^ pt > 20 GeV. 
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